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Abstract 

A real tunneling solution is an instanton for the Hartle-Hawking path integral with vanishing 
extrinsic curvature (vanishing "momentum") at the boundary. Since the final momentum is fixed, 
its conjugate cannot be specified freely; consequently, such an instanton will contribute to the 
wave function at only one or a few isolated spatial geometries. I show that these geometries are 
the extrema of the Hartle-Hawking wave function in the semiclassical approximation, and provide 
some evidence that with a suitable choice of time parameter, these extrema are the maxima of the 
wave function at a fixed time. 
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1. Introduction 



In the Hartle-Hawking approach to quantum cosmology the wave function of the 
universe is obtained from a Euclidean path integral over metrics and matter fields on a 
manifold M with a single boundary component dM. Such a path integral depends on the 
induced metric h a b and the matter configuration (p\gM on dM, thus determining a functional 

*[h,<p\ 8M ]= y £ l f[dg][M exp{-I E [g,ip]}, (1.1) 

M J 

where the summation represents a sum over topologies. The functional \& is the Hartle- 
Hawking wave function, which we are instructed to interpret as an amplitude for finding 
a universe characterized by h and <p\dM- The Hartle-Hawking approach neatly finesses the 
question of initial conditions by omitting any initial boundary, and it postpones the question 
of the nature of time in quantum gravity: information about time is hidden in the boundary 
geometry h, but the path integral can be formulated without making a choice of time explicit. 

In references || and it was shown that the sum over topologies in ( |1.1| ) can have a 
drastic effect on the behavior of the wave function. For most of this paper, this complication 
will be ignored; we shall focus on the contribution of a single fixed topology M to the wave 
function of the universe. In particular, we shall look for the maxima of the most probable 
spatial geometries and matter configurations of the universe. 

In the saddle point approximation, the wave function ( |1 . 1|) takes the standard form 

^M[h,(p\aM] ~ J2 D M[g,p]exp{-I E [g,Lp]} , (1.2) 

{g,<p} 

where the the sum is now over classical solutions {g, (p} with the prescribed boundary values, 
and the prefactor Dm is a combination of determinants coming from gauge-fixing and from 
small fluctuations around the extrema. In this approximation, it may not always be possible 
to specify h and (p\dM arbitrarily — for a given topology M, some boundary values may not 
lead to classical solutions — but we certainly have a great deal of flexibility. For instance, 
given one classical solution, we can cut off a small neighborhood of the boundary to obtain 
a new solution with the same topology but different boundary data. A given manifold M 
will therefore contribute to the wave function over a wide range of values of h and <p\dM- 

A particularly interesting class of saddle points consists of the so-called real tunneling 
geometries ||, solutions of the field equations for which the metric has Riemannian signature 
in M and vanishing extrinsic curvature at the boundary dM. Classically, these are the 
solutions for which the metric can be smoothly joined to a metric with Lorentzian signature 
to the future of <9M.[] Real tunneling geometries clearly give the right classical picture of the 
Hartle-Hawking "no boundary" boundary condition, describing universes characterized by a 
Lorentzian metric "now" but no initial boundary. Moreover, by analogy with more tractable 
problems such as pair production in a constant electric field , we expect these solutions to 
play an important role in the overall behavior of the path integral. 

* "Smoothly" here means "with finite action." 
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Note, however, that a real tunneling solution will typically contribute to *&[h, (p\dAf\ & t 
only a few isolated values of h ab - In the Hamiltonian formulation of general relativity, the 
extrinsic curvature is canonically conjugate to h ab , and we should not expect to be able to 
freely specify h ab while simultaneously maintaining K ab = 0. This conclusion can be checked 
explicitly in three dimensions with <p — and A < 0: a real tunneling solution is then a 
hyperbolic three-manifold with a totally geodesic boundary, and for any given topology M, 
the boundary metric h ab either does not exist (if there is no classical extremum) or is unique. 

We shall see below that the boundary values of real tunneling solutions are in fact the 
extrema of the wave function In general, these extrema are saddle points rather than 
maxima. But we shall find some evidence that with a suitable choice of time parameter, the 
real tunneling solutions determine the genuine maxima of the Hart le- Hawking wave function 
at a fixed time. 



2. Scalar Fields 

Before tackling the more difficult problem of gravity, it is instructive to look at the 
behavior of the "Hartle-Hawking wave function" for a free scalar field. The Euclidean action 
for a field ip in n spacetime dimensions is 

I E [<p) = ~ / d n x yfg (g ab d a pd b tp + m 2 p 2 ) , (2.1) 

Z JM v ' 

where for the moment the metric g is fixed. The first variation of the action is 

8I E [(p] = -f d n x^/g(5(p)(A g -m 2 )ip + j d^x Vh (8(p)V n <p, (2.2) 

J M JdM 

where V n is the normal derivative at dM and A g = V a V a is the Laplacian. Solutions of the 
classical equations of motion are thus genuine extrema as long as ip is fixed at the boundary. 



For a free scalar field, the path integral (1.1) is trivial; the saddle point approximation 



(|1.2|) is exact, with Dm = det(A 9 — m 2 ) -1 ' 2 . In this context, equation ( |2.2|) now has a new 
interpretation: it tells us that if we vary the boundary data <p\sm, the classical action IeW\ 
is extremal when V n ^ vanishes. This is the scalar analog of the condition for a real tunneling 
solution — a classical solution (p has a finite-action extension across dM to a Lorentzian 
spacetime precisely when V n <p = 0. We thus see from (|1.2|) that the scalar Hartle-Hawking 
wave function has its extrema at precisely the boundary values of the "scalar real tunneling 
solutions." 

We can now ask whether these extrema are maxima. If we restrict our attention to the 
space of classical solutions, the second variation of the Euclidean action, evaluated at the 
point V n ip = 0, is 

6 2 I E [<p]= I d n - 1 xVh(5^)V n (8p)= I d n x^(g ab d a (5^d b (8p) + (5^A g (5^)). (2.3) 

JdM JM v ' 

But for variations in the space of classical solutions, 

(A g -m 2 )(5p) =0, (2.4) 
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5 2 I E [<p] = / dfr y /g(g*d a (6<p)d b (6<p)+m 2 (6<p) 2 ) > 0. (2.5) 
Jm v ' 

Scalar real tunneling solutions are therefore minima of the Euclidean action, and hence 

maxima of the wave function. 

The scalar Hartle-Hawking wave function is thus peaked at the boundary value ip\dM f°r 

which the extension (p to M has a vanishing normal derivative. The map from (p\dM to V n <p 

is known as the Poisson map and we can summarize our results by stating that the 

maxima of the Hartle-Hawking wave function occur at the zeros of the Poisson map. For 

the more complicated case of gravity, the analog to (p\dM is t ne boundary metric h ab , while 

the normal derivative V n </? corresponds at least roughly to the extrinsic curvature K ab , and 

the Poisson map takes the metric boundary data to the corresponding extrinsic curvature. 

Let us now ask whether our results for the scalar field can be extended to this situation. 



3. Gravity 

For simplicity, we shall concentrate on the case of pure gravity, with a nonvanishing 
cosmological constant acting as a stand-in for matter fields. The Euclidean action I E is 

I E [g] = --^ I d n x sfg (R[g] - 2A) - ^— f dThcVhK, (3.1) 
16itG Jm 8tvG JdM 

where R[g] is the scalar curvature, A is the cosmological constant, and K is the trace of the 
intrinsic curvature of dM. Let n a be a unit vector field in a neighborhood of dM whose 
restriction to dM is the outward unit normal; for simplicity, we can choose n a to satisfy 

n b V b n a = (3.2) 

near dM. By a slight abuse of notation, let 

Kb = g a b - n a n b . (3.3) 

(The restriction of h a b to dM is the induced spatial metric on the boundary.) An easy 
extension of Wald's calculation in || then gives 

- lQnG5I E [g] = f d n xj-g {G ah + Ag ab )5g ab 

+ f d n ~\ \{K ab - h ab K)5g ab + h ab V a (n c 5g bc + 2g bc 5n c )] ■ (3.4) 

JdM L J 

Now consider the vector v b = n c 5g bc + 2g bc 5n c in a neighborhood of dM. Using the fact 
that n a is a unit vector, one may easily check that n a v a = 0, so v a is tangential on dM. Thus 

h ab V a v b = D a v a , (3.5) 



where D a denotes the covariant derivative on dM . The last term in fl3.4|) is therefore the 
integral of a total derivative on dM, and hence vanishes, giving 

~16nG5I E [g]= f d n x sfg (G ab + Ag ab )5g ab + f d n ~ l x Vh (K ab - h ab K)5g ab , (3.6) 
Jm JdM 
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the gravitational version of equation 

As in the scalar case, this relation tells us that solutions of the classical equations of 
motion are genuine extrema as long as the metric is fixed on dM. But as in the scalar 
case, it also tells us more: if we vary the boundary data, the extrema of the classical action 
are precisely the real tunneling solutions, for which K ab = on dM. For the scalar field, 
the saddle point approximation was exact, and these extrema were also the extrema of the 
Hartle-Hawking wave function. For gravity, this is no longer true — the prefactor Dm in 
( |1.2| ) can depend nontrivially on the boundary data, and there will be higher loop corrections 
- but as a first approximation, the real tunneling solutions will again be the extrema of 
*&[h], with corrections suppressed by powers of Planck's constant. 

To determine whether these extrema are maxima, we must again compute the second 
variation of the Euclidean action at the point K ab = 0: 

WirG5 2 I E [g\ = f d n ~ l x \^ 7 a6 5(# a6 - h ab K), (3.7) 

J dM 

where the variation of the metric has been denoted by 5g ab = lab (note that 5g ab = — r ) ab ). 
A set of rather routine calculations is described in the Appendix; the end result is that with 
the gauge choice [n a ^ ab h b c )\ dM = 0, 



32nG5 2 I E [g] = [ d n - l xVh 

J8M 



(3.8) 



where 7 = g abr f a b and (3 a = V b ( , y ab — \g ab l)- This relation is of the same general form as 



( |2.3j ), but in contrast to the scalar case, it is not manifestly positive. In particular, let us 
assume that we can consistently choose harmonic gauge, f3 a — (see the Appendix for a 
brief discussion of this choice). Then using Lichnerowicz's results for the variation of the 
Ricci tensor JTO]], 

5R ab = -]-A g j ab + ~ (~f a c R bc + ~f b c R ac ) - Y d Racbd = ^^7a6, (3.9) 



we obtain 



32nG5 2 I E [g] = J M dnx V9 



V c7 a b V c 7 a6 -^V c7 V c 7 



4A 

(n-l)(n-2) 



^7 7 Cacbd + — 7777 ^7 ( 7afe7 + —7 



(3.10) 



which has no clear positivity properties. 

Indeed, it is easy to see that the boundary value of a real tunneling solution is generically 
a saddle point of the Hartle-Hawking wave function, and not a maximum. For an empty 
space solution of the gravitational field equations with A ^ 0, the Euclidean action I E is 
proportional to the volume of M. But we can always decrease this volume by cutting off a 
small neighborhood of the boundary dM, or increase it by extending M slightly past dM. 

Intuitively, though, such a movement of the boundary is not the kind of deformation we 
should be concerned about. It is well known that the spatial geometry described by the 
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metric h a b implicitly contains information about time |TTJ, and a movement of the entire 



boundary dM corresponds physically to a time translation. It is not so surprising that the 
wave function has no maxima if such translations are allowed; the more interesting question 
is whether \l/ [h] has any maxima at a fixed time. 

To answer this question, we must extract information about time from the spatial metric 
hab an d the extrinsic curvature K a ^ at the boundary dM. This is not easy to do, since there 
is no unique way to parametrize time, but the structure of the real tunneling geometries 
suggests one natural approach: we can use York's "extrinsic time" K as a time parameter 
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If we restrict ourselves to "fixed time" variations of the metric, those for which K remains 
zero, we can eliminate at least some of the negative terms from (|3.10| ). From equations (|A.9|) 
and (|A.12 ) of the Appendix, we see that in harmonic gauge, 5K = implies that 



V n7 = 0. 

On the other hand, equation ( |3.9|) tells us that 

4A 

A ' 7+ (^2) 7 = - 



(3.11) 



(3.12) 



For a metric g a b with Riemannian signature, the Laplacian A g with boundary conditions 
( 3-1 1| ) has no positive eigenvalues. Hence for A < 0, 7 must vanish, eliminating the term 
V C 7V C 7 from (|3.10|) . For A > 0, equation ( |3.12| ) may sometimes have nontrivial solutions, 
but the spectrum of the Laplacian is discrete, so there will be at most finitely many of them; 
for most variations of the metric, 7 will again vanish. 

The term in ( |3.10|) involving the Weyl tensor is harder to control. To temporarily avoid 
this problem, it is interesting to look at the simple model of quantum gravity in three 
spacetime dimensions, where C a b c d is identically zero. Here, at least, it may be possible to 
prove that 5 2 Ie is strictly positive when SK = 0. 

Suppose first that A > 0. We must then check two conditions: that harmonic gauge 
P a = 0, assumed in the derivation of (|3.10| ), can always be reached, and that (|3.12| ) has no 
solutions with 7 7^ 0. If these conditions are satisfied, the second variation ( 3.10| ) will be 
manifestly positive. 

Both conditions require knowledge of the eigenvalues of Laplacians on M, acting on 
functions in (|3.12| ) and on one-forms in ( |A.6|) . Since these Laplacians are defined in terms 
of a background metric g a b for which dM is totally geodesic, they extend smoothly to the 
(closed) double M of M formed by attaching two copies of M along their boundaries. The 
boundary conditions ( [3. lip and ( |A.5| ) then guarantee that any solution of ( |3.12j ) or ( ]A.6p 
extends to a solution of the corresponding equation on M. Moreover, any closed three- 
dimensional Einstein space with A > is a quotient of the three-sphere by some group 
of isometries 
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so it is sufficient to look at Laplacians on the sphere. The appropriate 
eigenvalues are then given by Ray [14|, and it is not hard to check that both of our conditions 
are met. 

If A < 0, our task is somewhat harder, since the last term in (|3.10|) is no longer positive. 
In this case, the question can be reformulated in terms of hyperbolic geometry: the Euclidean 
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action is proportional to the volume of M , and the problem is to show that among hyperbolic 
metrics on M with K = on dM, the metric for which dM is totally geodesic has the smallest 
volume. One possible approach is the following. Hyperbolic structures on M — metrics with 
constant curvature —1 — are parametrized by a finite-dimensional moduli space M.. It is 
plausible that the volume of M becomes large outside a compact region of A4, where the 
metric on dM starts to degenerate. If this is true, and if the volume is a reasonably well- 
behaved function on A4, then it must have a minimum. But among hyperbolic metrics on 
a fixed manifold M, the real tunneling solution is unique, and is the only extremum of the 
volume; it must therefore be the minimum. It should be stressed that this argument is at 
best a strategy for a proof; work on this question is in progress. 

In four dimensions the situation is more complex, since the Weyl tensor in ([3.10 ) need not 
vanish. The inclusion of matter is likely to provide an added complication. For the moment, 
it seems difficult to make any definitive statements about the maxima of the Hartle-Hawking 
wave function for realistic cosmologies. But the partial results from three dimensions suggest 
a conjecture: that the maxima of the Hartle-Hawking wave function at "time" K — are 
precisely the boundary values of the real tunneling solutions. 



4. Conclusion 

The Hartle-Hawking program offers an intriguing approach to quantum cosmology, but 
it has not been easy to extract predictions from the formalism. Up to now, most work has 
involved either simple minisuperspace models (for example, |15[]) or classical or semiclassical 
approximations P,^ JTB| , ^7|] . It is therefore of some importance to understand how good these 



approximations are. 

If the boundary values of real tunneling solutions are the maxima of the Hartle-Hawking 
wave function, then classical and semiclassical approximations may be reasonably reliable. In 
particular, the classical signature-changing solutions of the Einstein equations are precisely 
the real tunneling geometries, and the interesting predictions based on these solutions [[17]] 
presuppose that they dominate the quantum wave function. Similarly, investigations of the 
sum over topologies [|,|],[16| have been based on the properties of real tunneling solutions, 
and may be invalid if the peaks of the wave function lie elsewhere. On the other hand, if the 
real tunneling contributions are not the maxima of the Hartle-Hawking wave function, this 
would have interesting implications as well: the wave function would then have no maxima, 
and would presumably be dominated by behavior near some boundary of the space of spatial 
geometries. 

This paper has not answered the question of whether the Hartle-Hawking wave function 
has maxima. We have seen, however, that if any such maxima exist, they must come from 
the contributions of real tunneling geometries. In the simple model of quantum gravity 
in three spacetime dimensions, we have found some reasonably strong evidence that these 
extrema are indeed maxima, but in realistic four-dimensional gravity the question remains 
open. 
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Appendix 

In this appendix, we fill in some of the details in the derivation of equation ( |3.8|) from 



(|3.7j ). Note first that the tangent directions to dM are defined without reference to the metric 
on M, so the direction of normal n a is also independent of the metric. We can therefore 
write | 

5n a = an a = -(^ bc n b n c )n a on dM. (A.l) 

To proceed further, it is useful to partially fix the gauge. An infinitesimal diffeomorphism 
corresponds to a transformation 

lab -> lab + V £ 6 + V b £ a , U^a = OU dM, (A.2) 

and metrics that differ by such transformations are physically equivalent. We can use this 
freedom to set 

n a j ab h b c = on dM, (A.3) 

implying that 

~f ab n b = 2an a on dM. (A.4) 
The remaining diffeomorphisms are now restricted to those that preserve ( |A.3| ), i.e., 

n a Ca = 0, V n (£^ a c ) = OI1 dM. (A.5) 

(We have used the fact that K ab = 0, which together with ( |3.2| ) implies that V a n b = on 
dM.) 

We could fix the remaining freedom by using the fact that 

Pa = V b (lab ~ habl) - Pa + (A„£a + R a %). (A.6) 



Given the boundary conditions (|A.3| ), V a b = A g 5 a b + R a b is a self-adjoint elliptic operator. 
With appropriate assumptions of smoothness, this means that the equation T> a b ^ b = — (3 a 
can be solved, thus transforming (3 a to zero, unless /3 a is itself a zero-mode of V a b . Such 
zero-modes are rare: in the case of empty space with a negative cosmo logical constant, for 
instance, the eigenvalues of T> are strictly negative. To postpone dealing with the exceptional 
cases when D has zero-modes, however, we shall not yet impose any condition on f3 a . 
We can now evaluate the variation of the extrinsic curvature K ab : 

5K ab = h a c VMb - \K c n d (V c i bd + V b i cd - V d i bc ) 

= -h a c V b (an c ) + \h a c Wni bc = ±h a c V n -f bc , (A. 7) 
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where the fact that V a n b = on dM has been used repeatedly. Thus 

^ ab 5K ab = |( 7 afe - n a n cl cb )V nlah = § 7 ab V n7o6 - 2aV„a, (A.8) 

and 

5K = g ab 5K ab = ±V n7 - V n a. (A.9) 

Note also that 

Kbl ab = 7 - n a n bl ab = 7 - 2a. (A.10) 
Inserting these results into equation ( |3.7|) , we find that 

32kG5 2 I e = [ d n ~ l x Vh~h ab V nla b - 7V n7 + 2 7 V n a + 2aV n7 - 8aV n al . (A. 11) 

JdM L J 

This equation can be further simplified by noting that 

n a (3 a = n a V b j ab - - V n7 = n a n b V nl ab + n a h b c V c j ab - ^V n7 

= 2 V n a + W c (2an b ) - - V n7 = 2 V„a - - V n7 , (A. 12) 

leading directly to equation (|3.8j). 
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